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Abstract 

The decay of Burgers turbulence with compactly supported Gaussian "white noise" 
initial conditions is studied in the limit of vanishing viscosity and large time. Probability 
distribution functions and moments for both velocities and velocity differences are com- 
puted exactly, together with the "time-like" structure functions T„(f, r) = {{u{t + t) — 

u{t)r). 

The analysis of the answers reveals both well known features of Burgers turbulence, 
such as the presence of dissipative anomaly, the extreme anomalous scaling of the velocity 
structure functions and self similarity of the statistics of the velocity field, and new 
features such as the extreme anomalous scaling of the "time-like" structure functions and 
the non-existence of a global inertial scale due to multiscaling of the Burgers velocity 
field. 

We also observe that all the results can be recovered using the one point probability 
distribution function of the shock strength and discuss the implications of this fact for 
Burgers turbulence in general. 



1 Introduction. 



The study of decaying Burgers turbulence (DBT) is largely motivated by the the observation 
that this is a system which falls into the phenomenological class of turbulent systems which 
can be treated in principle by means of Kolmogorov theory. Yet the answers which can be 
derived analytically for Burgers turbulence are in the sharp contradiction to the predictions 
of Kolmogorov theory. The understanding of the reasons for such a discrepancy and their 
relevance for the general theory of turbulence is one of the major aims of the study of Burgers 
turbulence. 

* tribe@maths.warwick.ac.uk 
^olegz@maths. Warwick. ac.uk 
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The history of the subject (see e.g. Q, [H, @, H, |6|, @, H], §, i, Q, 
|], H], H], 13], m, g, @, m, H, §; see for a review) shows however 



that the problem is hard, so hard in fact that it has a tendency to become self justifying, 
getting more and more alienated from the main body of turbulent research. However, until 
recently there existed no model of Burgers turbulence which can be used as a testing ground 
for general phenomenological theories of turbulence on one hand and admits a complete and 
simple analytical treatment on the other. 



In present paper we introduce and analyse such a model. Namely, we study the decay of 
Burgers turbulence with compactly supported Gaussian "white noise" initial conditions. In 
physical terms the turbulence in our model is excited by an initial disturbance localized at 
a fixed scale much less than the size of reservoir and which can occur with equal probability 
around any point of the reservoir. Note that DBT driven by "white noise" plays a special 
role for the theory of DBT in general. The reason is that integral scale of turbulence in this 
problem is not imposed by initial conditions but rather is generated by time evolution. Thus, 
the answers one obtains for "white noise" DBT are in some sense universal. Consider for 
example DBT driven by Gaussian initial conditions characterized by the two point function 
xir) which is approximately constant for r « R and goes to exponentially fast for r » 
R. Then the statistics of the velocity field in this model at scales much larger than R 
and much less than the integral scale is asymptotically equivalent, in the limit as ^ 
0, t — > oo, to that of "white noise" DBT. Likewise, compactly supported "white noise" DBT 
defines a universality class of models of DBT driven by compactly supported Gaussian initial 
conditions. 



The choice of a simple initial condition and the choice to look for answers only in the 
vanishing viscosity and large time limits lead to a model that is exactly solvable. Explicit 
asymptotics can be obtained for statistics that are hard to estimate in more general models. 
The main reason for the exact solvability of our model is the fact that the statistics of the 
velocity field in the case of compactly supported initial conditions are dominated in the limit 
^ 0, i ^ oo by two shock configurations, the statistics of which is easily computable as 
functionals of white noise. 



We would like to stress that our model is in a different universality class than the 
original Burgers model in which turbulence is initiated by white noise initial conditions but 
no restriction of compactness is imposed: a solution to Burgers equation corresponding to an 
initial condition supported on a whole line will generically contain infinitely many shocks at 
any moment of time, not just two as in our case. Accordingly, the large time statistics of the 
velocity field in our case is very different from that in Burgers' model. For instance, energy 
density decays as in our case (see section 3.1) and as t'"^/^ in Burgers', |11]. 



The paper is organized as follows. In section 2 we give a precise statement of the problem, 
construct a large time limit of the solution to the inviscid Burgers equation corresponding to 
compactly supported initial conditions and formulate the main statements about the statistics 
of these solutions. In section 3 we obtain asymptotics for a variety of statistics: the moments 
of velocity field; the probability distribution function of velocities; the velocity structure func- 
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tions; the probability distribution function of velocity differences; time- like velocity structure 
functions. In section 4 the analysis of these results is given. In particular, the validity of one 
shock approximation and multiscaling in the problem are discussed. 



2 The limiting velocity field. 



Consider the following initial value problem connected to the Burgers equation: 
du du d'^u 

— (x, t) + u{x, t)-Q^{^-, t) = ^'q^^^-' ^ ^ * > 0' (1) 

n(x,0) = uq{x), (2) 

where uq{x) is a bounded function which is compactly supported in the interval [xq — 1, xq + 1]. 
Here / is a fixed positive constant and xq is a random variable uniformly distributed in 
the interval [— L,L]. The fixed positive constant L plays a role of normalization length. 
Conditional on xq the initial velocity uq{x) will be a white noise over the interval [xq — xq+Z], 
so that it has a formal density 

^'(no|xo) = |e'^C''/"o(")'^ (3) 

where Z is a normalization constant chosen in such a way that, formally, 

'■^ dxQ 



j P{n,\x,)D{u,) = l. 



J, the Gaussian variance, is a positive constant which plays a role of Loitsansky integral for 
the problem at hand. 

Since we have a compact initial condition the distribution of the velocities ^((x) are not 
translation invariant. The role of xq is to randomise the location of the initial disturbance 
uniformly over the interval [— L,L]. The values of ut{x) at a fixed x will then typically be 
non-zero only with probability 0{L~^). We take the limit as L — > oo and all the answers 
concerning the statistics of the velocity field will be expressed in the form of the leading 
term in an asymptotic expansion in L^^ . This has the advantage that the answers are then 
translation invariant and we are free to consider statistics centered at the origin. 

In what follows we will compute asymptotics of the following statistics: the moments of 
velocity distribution M„ = (u"(x,t)); the velocity structure functions Sniy) = {{u{x + y,t) — 
u{x, i))"), the probability distribution function of the velocity field P{u) = {6{u—u{x, t))); the 
probability distribution function of velocity differences P{u,y) = {6{u — u{x + y^t) + u{x,t))); 
and the "time-like" velocity structure functions, Tn{T,t) = {{u{x,t + r) — M(x,t))"). Here 
6{z) = x{z ^ 0) is the Heavyside function and (. . .) denotes the average w.r.t. to the random 
initial velocity field uq[x). 

The solution of the initial value problem (||), (|2|) for i/ > via the Cole-Hopf transfor- 
mation and the evaluation of the limit as — > for fixed t > are well known. We refer the 
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reader to |20f| and 1 11 ] for a detailed description and give here a quick summary, sufficient for 
our needs. The vanishing viscosity solution can be obtained by plotting a chain of parabolic 
arcs such that each is touching the graph of the function —q{x) = — UQ{y)dy at two points 

exactly. The i-th parabolic arc is given by a graph of the function ^i{x, t) = + As 
time grows the parabolic arcs flatten out and merge, and there exists a time T* such that for 
any t > T* there are generically only two arcs left. The velocity field associated with such a 
configuration is then given by 

u*{x,t) = U{xo + x%x,t,P,Q) = ^ X[(,„+,._^/32^,,„+,«+^/5(;^ZQj7](a^), (4) 

where Xi is an indicator function of the interval I, P = q{+oo) is a momentum corresponding 
to a given uq, Q = minx q{x) is a global minimum of q{x) and xq-\- x* G [xq — /, xq + I] is the 
point where this minimum is achieved. (Such a point exists and is unique almost surely as 
q{x) is continuous and the global minimum is almost surely unique.) The limiting solution 
(0) was originally constructed in |20|. 



The time T* at which the limiting velocity field u* is attained depends on the random 
initial condition uq but it will be shown that the statistics of the velocity field is well ap- 
proximated at large times by the statistics of the limiting velocity field u* . The latter is 
determined in turn by the joint distribution of the momentum P and the global minimum Q. 
Indeed although the expression for u* depends explicitly on (P, Q, x*), the dependence on x* 
doesn't influence the statistics of u* in the limit L — > oo, where the translational invariance 
is restored. We delegate the detailed discussion of this point to the next section. 



The choice of white noise as an initial distribution leads to the distribution of the pair 
(P, Q) being exactly calculable. Indeed it is a well known consequence of the 'reflection 



principle' for Brownian paths |29]. Since it is key to all our asymptotics we include a quick 
derivation of the joint density function p{P, Q). We start with a computation of the proba- 
bility distribution function of momentum p{P). Writing 5 for the delta function at zero, we 
have by definition 



p{P) = {5[P- I dxuo{x)])= I ^—e'^^{e 



Using (^) this functional integral is Gaussian and can be simply computed to give 



e 

The integral over A is a Gaussian integral and we conclude that the distribution of P is also 
Gaussian, as could have been guessed from the very beginning, and given by 

p{P) = ^7^, where Pq = ^VU . (5) 

The joint probability distribution function can now be computed as follows. Fix q,p satisfying 
q < 0,q < p. Let x' be the first value of x for which q{x) = q. Define q'{x) to equal q{x) 
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for X < x' and to equal the reflection of q{x) in the horizontal line y = q for x > x' . Then if 
Q' = m.inxq'{x) and P' = q'{oo) the reflection principle (see p9|), which exploits the white 
noise nature of uq, states that Q',P' have the same distribution as Q,P. Then 



Prob(Q <q,P>p) 



Prob(Q' <q,P' <2q- p) 
Prob(P' <2q-p) 




Differentiating in p and q we conclude that 



4(p - 2g) -(£p22)2 



if g < ■min{0,p} 



(6) 



and is zero for all other values of p and q. 

With the help of (P) we are able to average functional F[u*{t)] = F[u*{xi,t) : i = 
1,2,.. .] with respect to the initial distribution. If however we are interested in the statistics 
of u{x, t) at zero viscosity and large times there is still a question: is it true that in this limit 
{F[u{t)]) {F[u* (t)]) , or even at large times are there statistically many initials conditions 
such that corresponding velocity profiles haven't converged to the limiting ones? It so happens 
that the first alternative prevails. The detailed proofs of this fact for relevant functionals are 
carried out in the next section and in the appendix and are based on the following estimate 
on the time T* of convergence to the limiting profile: 



where ic = y -j and C is a positive number. The proof of this estimate is fairly complicated 
and is allocated to the appendix. However the result itself is so important for the validity 
of conclusions of our paper that we decided to present here a convincing and very simple 
heuristic derivation of it. 

By definition, Prob(r* < t \ P,Q,x*) = Prob(g < \ P,Q,x*), where coincides for 
X < x* with with parabolic arc passing through the point {x*,—Q) and touching the 
line y = and with parabolic arc ^i^t passing through the point (x*, —Q) and touching the 
line y = —P for x > x* (we used the translation invariance of the random variable T* to set 
xq = 0. Consequently, x* £ [—1,1]). 

It is convenient to think of a Brownian walk q(x) passing through (x* , —Q) as a collection 
of two independent walks q~^{x) and q~{x) starting at this point and moving in the opposite 
directions in "time" x. Therefore, 

Prob(T* < t I P, Q, x*) = Prob(g- < | Q, x*) ■ Prob(g+ < $2,t I Q, x*). (8) 

To estimate, say, Prob(g^ < ^i^t \ Qi^*) below we note that Prob(g'^ < ^i^t \ Q,x*) > 
Prob(g'^ < —Q + 9 ■ {x — x*) \ Q,x*), where y = —Q + 6 ■ {x — x*) is an equation for the line 




(7) 
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tangent to the parabola at the point (x*, —Q); 



Hence, 



Prob(g- <^i,t\ Q,x*) > 



> hm 



q{x*)=-Q-e 

J DqQ 

q(-l)=0 



q < -Q + 6 ■ {x - X*) 



q{x*) = ~Q-e 
q(-l)=0 



q<-Q 



e O<-Q-0-(/ + x*) , (9) 



where 0[. . .] is a functional step function, 0(. . .) - a usual one. 

The functional integral in the numerator of (|9|) can be transformed into an integral over 
all pathes satisfying q{x) > by a change of variables q{x) — > q{x) — Q + 9{x — x*). (A counter 
part of this transformation in quantum mechanics is a Galilean transformation.) Now the 
functional integrals in both numerator and denumerator of (P) can be expressed in terms 
of Green's function of heat equation q = ^q" on half a line, i.e. the antisymmetrization of 
Green's function of the same equation on the whole line. A simple computation shows then 
that 



Prob(Q- < $i,t \Q,x*)> 1 



2/2 



-Qt 



2f 



-Qt 



(10) 



Similar estimate holds for Pioh(q~^ < $2,* I P, Q, x*) if one replaces —Q with P — Q in the r. 
h. s. of (lie 



Substituting these two estimates into (^) and integrating both sides of the resulting 
inequality w. r. t. P, Q using (^) we find that Prob(T* < t) > 1 — Const^ which is 
equivalent to the estimate (0) for Prob(r* > t) = 1 - Prob(T* < t). 



3 The statistics of the velocity field in the 1/ ^ 0, t ^ 00 limit. 
3.1 Moments of the velocity distribution. 

The aim of the present section is to compute the large t-limit of moments of the velocity 
distribution 



M„(t) = (n"(0,t)), n = 1,2,... 



(11) 



Odd order moments vanish identically due to the symmetry: both Burgers equation and the 
initial distribution are invariant with respect to the transformation u —u, x — > —x. On the 
other hand, M2k+i — > —M2k+i under this transformation, which implies that M2fc+i(t) = 
for fc = 1, 2, . . . We concentrate therefore on the computation of the moments of even order and 
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assume everywhere below that n is even. We may write, using the fact that u{x, t) = u*{x, t) 
for t > T*, 

Mn{t) = {u*''{0,t))+Rnit), (12) 

where 

Rnit) = ^ (^^."(0, t) - n*"(0, t)^ e{T* - t)J (13) 
is an error term to be estimated. 



The first term in the right hand side of (12) can be written in the following form: 

'■'^ dxQ 



{u*^{0,t))= J dpdqp{p,q) |^^^C/"(xo,0,t,p,g)+r„(t), (14) 



where 



rn{t) = dx* J dpdq p{p, q, x*) J_^^^^ + J^^^'^ ) ^ [/'^(xq, 0, t,p, q) (15) 

is an error term appearing due to neglecting x* in comparison to L and p{p, q, x*) is a joint 
probability density of P, Q and x* . It is shown in the appendix that the error term r„(f) does 
not affect the asymptotics as ^ 0, L ^ oo, t — > oo. Informally this fact can be explained by 
noticing that the integrand in (|l5|) is non-zero only for velocity profiles which are "stretched" 
over the interval of length L and thus are exponentially improbable. 

The remaining integral on the right hand side of ( p!^ ) can be evaluated exactly using the 
explicit expressions @) and (^) leading to the following result: 

Vvr(n + 1) L 

where 

L{t) = Vm^t, U{t) = ^ (17) 

are parameters, with dimensions length and velocity, which should be interpreted as the scale 
of turbulence and turbulent velocity correspondingly. Here we write the symbol ~ to mean 
asymptotic equivalence in the limit as L — > oo and then t — > oo. 

Another computation presented in the appendix leads to the following estimate of the 



error term Rn{t) from (12): 

1/4 



Rn{t)\<Cn^U{tr(^^^ , (18) 



where tc = \/l^JJ is a constant having a dimension of time, C„ is a positive constant. 
Comparing ( [T^ ) with (|T6|) we see that for t » tc, \Rn(t)\ « (n*"'(0,t)), which permits us 
to conclude that 

r(fc/2 + 3/4) L{t) 

^''^'^^ V^{2k + 1) -T^^') ' ^ = 1'2,... (19) 
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It is important to stress however that coefficient C„ from (|18| ) grows faster with n than 
the number factor in the r. h. s. of ( |T9[). Thus it takes a long time for a moment of high 
order to converge to the limiting value (|l9|) . 

It follows from ( p^ ) that the energy density E{t) = ^M2{t) decays like as t ^ oo. 

This is the result to be expected: Dissipation of energy occurs in Burgers turbulence due to 
shock collisions and at each separate shock. The energy of a separate shock decays as 
and due to the absence of shock collisions in the limiting profile (^), this also gives the law 



of decay of total energy density. This argument is due to J. M. Burgers, see jll]. 



We will also see below that the statistics of the velocity field in our model is self-similar 
with the scales of length and velocity given by (|l^. These scales depend on time exactly 
as their counterparts in Kida's model. The statistics of the velocity field in our case are 
however different from that of KidaQ. Thus we conclude that the self-similarity alone does 
not determine the large time asymptotics of the statistics of velocity field in DBT. Note also 
that E{t) decays in time, showing the presence of a dissipation anomaly in the model: the 
rate of energy dissipation does not vanish but converges to a finite non-zero limit when the 
viscosity approaches zero. 



3.2 The probability distribution function of velocities. 

In this section we will concern ourselves with computing the probability distribution function 
(PDF) of velocities given by 

P(u, t) = Prob(u(0, t)>u) = {e{u{0, t) - u)) . (20) 

Reasoning exactly as in the previous section we find that 

P{u,t) = {9{u*{0,t) -u)) +R{u,t), (21) 

where 

R{u,t) = (^(^9(^u{0,t) -u^ -e(^u*{0,t) -u^^ e{t* -t)'^ (22) 
is an error due to the replacement u ^ u*; 

{9{u*{0,t)-u)) = 

= 0{-u)+ J dpdqp{p,q) J^^^(^e{U{xo,0,t,P,q) - u) - e{-u)^ +riu,t), (23) 
where r{u,t) is an error due neglecting x* in comparison with L: 

r{u,t) = 

= J^^dx* J dpdqpip,q,x*)yj^^^+ ^^i9iU{xo,0,t,p,q)-u)-e{-u)). (24) 

^There exists no complete solution of Kida's model. Yet the answers which can be obtained within Kida's 
model are different from their counterparts in our model. 
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The reason that the term 9{—u) is added and subtracted is that, due to the averaging of the 
position of the initial condition over the block [— L,L], the velocity is typically zero and so 
the PDF is an 0{L~^) perturbation to 6{—u). 



An estimate of r(n, t) similar to that of the term r„(f) in section (3.1) shows that r(n, t) 
does not affect the final asymptotics. An exact calculation using the known density p(p, q) 
for the other terms on the right hand side of (p3|) leads to 



{e{u*{<d,t)-u)) ^e{-u) + 



m 

L 



da _„2/ \ 
J«2 V / 



where u = u/U{t). A computation performed in the appendix shows that 

\R{n,t)\<C^Cj)^/' 



(25) 



(26) 



where C is a positive constant. Comparing (26) with (25) we see that for t >> tc we 
have {9{u*{0,t) — u)) >> \R{u,t)\, with the last inequality being pointwise in u rather than 
uniform. We conclude that 



P{Uit)u,t)r^9i-u) + ^ r ^e-"'" ( - \u\]sgniu) 
If in particular ^ oo this simplifies to 

P{U{t)u,t) ~ d{-u) + J-sgn(u)^^. 



(27) 



(28) 



Note that the answer (p7|) for P(n, t) is self-similar with U (t) playing the role of the 
integral velocity scale. Note also that the form of P{u, t) is not Gaussian. This confirms the 
non-triviality of our model: the output (the strongly non-Gaussian statistics of the velocity 
field in the limit of small viscosity and large time) is not the same as the input (a trivial 
Gaussian distribution of the initial velocity field). This non-triviality will be re-emphasised 
in the consequent sections where it will be shown that the limiting statistics of the velocity 
field is intermittent. 



Finally we would like to make the following technical comment. Of course, the moments 
of the distribution (|27| ) are exactly those given by ([l9|). We could therefore try to compute 
the distribution (20) first and then argue that the moments of this asymptotic distribution 
coincide with the asymptotics of the moments of the actual distribution. Unfortunately the 
analysis of error terms within this approach becomes very involved. For this reason we have 
two separate computations, the asymptotics of the moments of the velocity distribution and 
the asymptotics of the velocity distribution itself. 
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3.3 Velocity structure functions. 



Now we will turn to the two-point statistics of the velocity field and compute asymptotics 
for the velocity structure functions given by 

Sniy,t) = (^(^u{y,t)-u{0,t)^ ^ n = l,2,... (29) 

We find as in the previous subsections that 

Sn{y,t) = (^(^u*{y,t)-u*{0,t)y^+Rn{y,t), (30) 

where Rn{y,t) accounts for the error due to the replacement of u with u* . As shown in the 
appendix this error can be estimated as follows: for t such that L{t) > y 

\Rn{y,t)\<CnU^{t)l(^^y\ (31) 

We express the first term in (^0|) as 

(^(u*{y,t)-u*{0,t)^ 

where r^iy, t) accounts for an error arising due to neglecting x* in comparison with L. Again 
it can be shown that the term rn{y,t) does not contribute to the asymptotics. Now a direct 
computation using the density p{p, q) shows, for y > 0, that 

u*{y,t)-u*iO,t)y) ~ (-l)-^T(^^]^U^t)y + Oif), n = 2,3,... 



vr V 4 7 L 



where y = j^- In addition Si{y,t) ~ 0, which confirms the restoration of translation 
invariance in the large L limit. Comparing this with (^) we see that the asymptotics of the 
velocity structure functions is given, for fixed y < 1, by 

Snimy, t) ~ (-l)"^r(^)^^7«(t)y + Oif) n = 2, 3, . . . (32) 



It has been assumed in our computations that y > 0. Extending (32) to negative y by the 
symmetry y — > —y, u ^ u, we see that S2k{L{t)y,t) is proportional to \y\ and S2k+i{L{t)y,t) 
is proportional to y for k >1 and \y\ « 1. 

Thus the velocity structure functions of the problem exhibit in the inertial range the ex- 
treme anomalous (non-Kolmogorov) scaling which is typical for Burgers turbulence in general 
and is due to the presence of shocks in the limiting velocity profile. The Burgers anomalous 
scaling is well known from heuristic arguments (see e.g. [|^], 0, 1^] ). In our case however 
it has been derived as a part of the complete solution of the problem. 



10 



3.4 The probability distribution function of velocity differences. 



Here we will compute the PDF for velocity differences 

P(n, y, t) = Prob (u > An(y, t)] =(e(u- Au{y, t)]), (33) 



where Au{y,t) = u{y/2,t) — u{—y/2,t) and y > 0. Definition ( p^ ) is tailored for the study 
of negative velocity differences and we consider only the case u < 0. Negative differences are 
the interesting case since they occur when the velocities are evaluated either side of a shock. 
A lengthy but straightforward computation shows that for fixed u < 0, y > 

P{U{t)u,L{t)y,t) ~ 
L(t) / fiy-"^)^ da ^2/ ^ \ /■°° da „2\ 



where, as shown in the appendix, 



\R{u,m<C^^^(Ty^'. (35) 



L y + \u\ \ t 

Due to the presence of extra factor of {^)^^^ decaying with time, R{u,y,t) becomes small 
compared to the first term in the right hand side of (p4[), given that u,y fixed. 



It is easy to analyse ( |34D in the following limiting cases. We suppose that y << 1. If 
|n| << 1, then 

P{U{t)u, L{t)y, t)^^[l-^y-^\u\^ + 0{f) + 0(|np) j . (36) 

If 1 << |n| << y^^/^ then 

P{U{t)u, Lit)y, t) ~ M^^e-I^l' f 1 + O(^) + Oiy\uf)) (37) 
If 1^1 >> y^^l'^ then 



To summarize, for negative P{u,y,t) decays algebraically for |n| << U{t) and super 
exponentially for \u\ >> U{t). Moreover, P{u,y,t) ~ 0{y) if 1 << \u\ << y~^^^ and doesn't 
depend on y if 1 << |n| << y~^^^. This information alone enables one to conclude that 
velocity structure functions of sufficiently high order exhibit anomalous scaling. In addition 
we observe a crossover between regimes (37) and (^). This crossover is actually responsible 
for the presence of many scales in description of the statistics of velocity field and the absence 
of the universal inertial range in Burgers turbulence. We refer reader to section 4 for a detailed 
discussion of this point. 
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3.5 The multi-time statistics of the velocity field. 



The simphcity of our model allows us to compute the correlation between values of the 
velocity field at different moments of time. Let 

TniT,t) = (^(u{0,t + T)-u{0,t)^ ^ (39) 

be the velocity structure functions corresponding to the same point at space but different 
moments of time. We write (|39|) in the already familiar form 



TniT,t) = (^(u*iO,t + T)-u*{0,t)^ ^+i?„(r,t) (40) 

with Rn{T) accounting for an error due to the replacement of u with u* . An estimate in the 
appendix shows that 



r 



U{t) ftc 



1/4 



\Rn{r,t)\<CnU-{t)^i^jj . (41) 

The computation of the first term in the right hand side of ( ^0[ ) is very close to the compu- 
tation performed in previous sections and leads to 

Y(ii±3\ TT(f)r 
T^{T,t) ~ (-l)«^^[/"(t)^ + i?„(r,i) 

U(f)T 

~ (-l)"4^[/"(t)^, n = 2,3,... (42) 

We therefore conclude that the time-like structure functions exhibit in Burgers turbulence 
the extreme anomalous scaling in r given by T„(r, t) ~ r, n = 2, 3, . . .. Comparing this with 
the expression (]3^) for the space-like structure functions, we see that 

Sn{y,t)=T^{T,t), n = 2,3,... (43) 



at y = C{n)U{t)T, given that y << L{t) and t «t. The identity ( pSj ) means that "isotropic" 
Taylor conjecture stating the equivalence of the space-like and time-like statistics in isotropic 
turbulence at small scales, becomes a theorem for our model of Burgers turbulence. The 
similar observation was also independently made in in the context Burgers turbulence 
generated by correlated Gaussian initial conditions. 

Let us finally note that if one wishes to compare T„(y, t) with Sn{T, t) at arbitrarily high 
orders n, the condition of applicability of relation ( |4^ ) has to be changed to y << L„(t), 
where Ln(t) is correlation length associated with n-th order structure function introduced in 
section 4.2. For n » 1, L„(t) ~ L{t)/n^/^, see M below. 



3.6 One-shock approximation. 



We wish to show that all of the results obtained in the previous section can be easily obtained 
from heuristic arguments given the knowledge of the probability density of a velocity jump 
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at a shock. In our case the latter is easy to compute: a simple computation which uses the 
knowledge of the limiting velocity profile (Q) and the density p{p, q) gives 



p[^)^(^5^^-^2{P-Q)/t^^ = -j=^^Jie-^^ , (44) 

where is a velocity jump at the (right) shock, Jl = -p^. The probability density of the 
velocity jump at the left shock has exactly the same form, so we will be referring to (44) as 
the probability density of the velocity jump at a shock. 

Now let us assume: Firstly that the large-f statistics of u are approximated by that 
of u*; secondly that a one-shock approximation is valid, i.e. that one can disregard in the 
analysis the contributions coming from configurations with shocks separated by distances 
much less than the average separation L{t). 

To derive P{u,y,t) for u < 0, y « L{t) using these assumptions note that u{y,t) — 
u{0, t) can be negative only if there is a shock at some point in [0, y]. If the right hand shock 
lies at X € [0,y] then u{y,t) — u{0,t) = —p + x/t. A similar formula holds if the left hand 
shock lies in [0,y]. So neglecting the contribution from the configurations with 2 shocks inside 
the interval [0,y], we see that 

/ \ py dx X 

Prob u{y, t) — n(0, t) < u] ^ 2 — Prob(Size of Jump > u). 

\ J Jo 2L t 

This can be easily computed using the density of the shock jump (Q) giving 

ProbUy,* -nO,t <2/U— ^ / ^e"'^ + n + ^ / 

\ J L \Ju2 yj-n: J{y-u)'^ ' 



which coincides with the exact answer 

With the knowledge of the PDF of velocity difference we can compute velocity structure 
functions, thus moments of velocities, thus the PDF of velocities. In other words all of the 
results of the previous section concerning single time statistics of the velocity field can be 
obtained using a one-shock approximation. 

Moreover, the r-dependence of the time-like structure functions (|^) is also entirely due 
to the one-shock effects: if n = 2,3, ... and r << t, then the main contribution to T„(r, t) 
comes from the configurations with a shock passing through x = between the moments of 
time t and t -|- r. A shock with velocity jump p travels a distance approximately [it j2 over 
the interval [t, i + r] . Therefore, 

r„(r) ^(— /i)"x^Shock passed through during [t, t -|- r] 



Computing this average using the PDF of shock strength (44) we arrive exactly at ([4^, which 



again shows that one-shock approximation is asymptotically exact. 
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These calculations support the following statement about decaying Burgers turbulence: 
all one needs to know in order to describe the statistics of the velocity field at scales much 
less than the average distance between shocks is the one-point PDF of shock velocity and 
strength (or just shock strength if the correlation functions which we're trying to compute 
are Galilean- invariant). Thus the problem is much simpler than one might have thought: 
recall for example that exact formulae expressing velocity correlation functions in terms of 
the statistics of shocks are such (|pA[|, p3[) that one seemingly needs to know the n-point 
joint PDF of shock strengths in order to compute n-th order correlation function. 



The rigorous proof of the above statement together with estimates on the errors of one- 
shock approximation will make DBT analytically tractable for a wide class of initial conditions 
as the great deal is known about the one-point function of shock strength, see e. g. 

m, a, a. 



Is there a universal technique for the computation of the one-point PDF of the shock 
strength? It has been known since Burgers [11|, but never really exploited, that shocks behave 
(almost) as a system of sticky particles. One might try therefore to extract the information 
about one-point PDF of shock characteristics by studying the kinetics of this system, for 
example, by analyzing the Smoluchowski-Bogoluibov chain of equations for one-point, two- 
point, . . . PDF's of shocks. 



3.7 On multiscaling in Burgers turbulence. 



In statistical physics the term "multiscaling", instead of "anomalous scaling", is used to 
stress an inherently multiscale nature of a system exhibiting anomalous scaling of correlation 
functions. Burgers turbulence is no exception. In this section we will show that the crossover 
between the tails ( |37| ) and (|3^) of the PDF for velocity differences is actually a reflection 
of the presence of many correlation lengths in the problem, which in turn is a consequence 
of the anomalous scaling of correlation functions and, ultimately, the intermittency of the 
velocity field in Burgers turbulence. 

Let n >> 1 be a large even positive integer. We know from (|3^ ) that as y approaches 

zero 

5„,(L(t)y,t)«-^r(!^)^C/«(t)y. 

For large y however one expects the quantities u{L{t)y, t) and u{0, t) to become independent. 
When this happens we have 

Sn{my,t) = ((n(L(t)2/,t)-n(0,t))") 
~ (n-(L(i)y,t)) + (n"(0,t)) 

^7r(n + 1) L 

Here the cross terms in expanding the nth power are, using the independence, of order 0{L~^). 
The region in between these two formulae for large and small y marks the correlation length 
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for the nth moments. If we assume there is a simple crossover then we can locate the scale at 
which it occurs by equating the expressions for large y and small y. These become equal, i.e. 
Sn{L{t)y,t) ~ 2M„, at the value n~^/^ and so the correlation length for the n-th structure 
function is 

^n~^, n»\ 46 
and this shows the presence of many scales in our problem. 

To show how this multiscaling is related to the crossover between the asymptotic regimes 



(pTl) and (35) we shall use the PDF for velocity differences to compute Sn{y,t) for n positive 
and large. Writing Sn{L{t)y,t) as an integral against the PDF of Au{L{t)y,t) and treating 
n as a large parameter we see that the integral is dominated by values of |u| coming from the 
neighbourhood of the negative critical point of the function 

F{u) = |urexp(-|u|^) 

namely near Uc = — n^/^. Note this value is much less than —1 for n >> 1 and so we 
may neglect the part of the integral that uses the PDF in the form (^) and also neglect 
positive values of Au{y). Now, if in addition \uc\ « y~^^^, we have to use asymptotics 
(|37|) to evaluate the contribution from the critical point, which yields Sn{L{t)y,t) ^ Cy. If 
>> y~^^^ we have to use asymptotics (38) in our computations, which gives Sn{L{t)y, t) « 



u, 



Constant. The crossover between these two answers corresponds to the crossover between the 
asymptotics (^6|) and (|37|) and occurs when y = \uc\^ = n~^/'^, exactly as in our computed 
correlation length L„ for the n-th structure function. 

It remains to remark that multiscaling, and consequently a PDF for velocity differences 
which has a crossover between a regime scaling like y and one that is independent of y, 
should be a general feature of DBT regardless of the initial distribution. All related questions 
concerning other statistics can be studied in more general situations, if one assumes a one- 
shock approximation is valid, by using the information about the tails of the one-point PDF 
of shock strength obtained in [32|, [^]. 



It is worth noting that the presence of the multitude of correlation lengths in Burgers 



turbulence was understood long ago by Robert Kraichnan, |25], and rediscovered within the 



instanton approach to the forced Burgers turbulence, |12|. It is also worth stressing that 
in models of chaotic systems which do not account for the effects of intermittency, there is 
always a single universal correlation length. A good example is served by random matrix 



models, see [26| for a review. 



Finally, let us remark that if we define the integral scale as the scale of scaling behaviour 
of correlation functions, we must immediately conclude that there is no such unique scale, 
there is rather a family of them parameterized by the order of correlation function. In other 
words, the notion of the integral scale becomes local, and the notion of the universal inertial 



range disappears. (See also [14| for the general discussion about the multitude of dissipative 
scales based on a multifractal models). This should be a general feature of all intermittent 
turbulent systems, for instance, Navier-Stokes turbulence. 
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5 Note added in proof. 



We are grateful to the referee of our paper who drew our attention to a recent preprint by L. 
Frachebourg and Ph. A. Martin, |13], in which the study of the model of decaying Burgers 
turbulence initiated by white noise initial conditions (without compactness assumption) has 
been effectively completed. This model was originally considered by Burgers himself about 
forty years ago but complexity of analysis prevented him from obtaining explicit answers 
for anything but the two- and three-point correlation functions of velocity field. Now most 
of the questions about the statistics of velocity field in Burgers' model can be effectively 
resolved using the integral representation of the Green's function of a diffusion equation in 
the {x, t)-domain with parabolic boundary derived in the above mentioned paper. 



6 Appendix. 



In order to bound the various error terms in section 3 we will need to bound the size of the 
true solution u, the asymptotic solution u* and the size of their supports (i.e. the interval on 
which they are non-zero). We use details from the method of construction of the vanishing 
viscosity solution as descibed in |20| and recalled in section 2. 



Suppose that initial velocity profile is supported in the interval (xq — /,xo + /)• The 
rightmost (respectively leftmost) parabola in the chain of parabolic arcs built on the initial 
potential will always lie to the left (respectively right) of the parabola with the same cur- 
vature that passes through the point (xq + /, —Q) (respectively (xq — /, —Q)) and assumes 
minimial value equal to —P (respectively 0). This immediately implies that both u and u* 
are supported in the interval y*] where 



= xo - Z - V^2tQ, y* =xo + l + ^j2t{P-Q). (47) 

Using the fact that both u and u* vanish at the point within [xq — /,Xo + I] at which q{x) 

^Who asked the very useful, perhaps rhetorical, question "Why study white noise Burgers turbulence at 
all?" 



16 



achieves its global minimum, we also find that \u\ and \u*\ are bounded by Umax where 

^imax = max{(y* - (xo - ((xo + /) - 

= max{{2l + ^f2t{p'^))/t,{2l + ^^2tQ)/t}. (48) 

Estimates (^7|), (^) and the bound (|^ will be used to estimate all relevant error terms in 
section 3. The careful analysis of these error terms leads to a better understanding for when 
the asymptotics for various statistics start to hold. 



6.1 Proof of the estimate ([TSl) . 

Applying the above estimates to the error term (^8|) we obtain 

\Rn{t)\ = (|n"(O,t)-n*"(O,t)|0(r*-t)) 

= (|n"(0,t)-n™(0,t)|x[,.,,.](0)e(T*-t)) 

< 2(uSiaxX[,.,r](0)^(^*-*)) 
2 

< -j-iiy* - y*)u^a^S{T* - t)) (averaging over xq) 

< jiiy* - y*)'<ax)'/'(^(r* - 0)'/' (Cauchy-Schwartz) 

< c„^[/"(^)(|)V^ 

where the last inequality uses the the estimate ^ and an explicit calculation using p{p,q). 
Comparing the first and the last entries of the presented chain of inequalities we obtain a 
proof of (18). 



6.2 Proof of the estimate on r„(t) from section ( |3.1|) . 

We may bound r„,(t) as follows: 

rl r r-L+x* rL+x* 

\rn{t)\ < J_dx* J dpdqp{p,q,x*){J^^ + )-l\u{xo,0,t,p,q)r 
< f dpdqp{p,q){r'^^\ f'^^')^\U{xo,0,t,p,q)r 
- Tvir) J dpdqp{p,q)ei,/^t-{L-l))U^{t) 

^ !m"(a^-(-(w)>"<" 

where the last inequality follows by an explicit calculation using p{p, q). This is exponentially 
small in L and so does not affect the asymptotics which take the limit L ^ oo first and 
preserve only the 0{L~^) terms. A similar argument controls similar error terms of this form 
for the other statistics considered. 



17 



6.3 Proof of the estimate ([261) . 



The proof of ( p6| ) is similar to that of ([T^): 

\R{u,t)\ < ={\e{u{o,t)-u)-e{u*{o,t)-u)\x[y„y^]{o)e{T* -t)) 

< 2(x[,„,.](O)0(T*-t)) 

< \{{y*-y*)Y'\0{T*-t))^i^ 

L ^ 



6.4 Proof of the estimate (|31]). 



We can spht this error term into two via 

\Rn{y,t)\ < 2^{Au{y,tre{T* - t)) + 2^{Au*{y,treiT* - t)), (49) 

where Au{y) = u{y/2,t) — u{—y/2,t). We show how to bound the first of these terms, the 
other being entirely similar. The vanishing viscosity solution u takes the form, within its 
support, of a line with slope 1/t plus a series of downward jumps. So we may define F{x,t) 
to be a non increasing piecewise constant function so that, for x in the support of 

u{y, t) = ^—p^ + F{y - xo,t). 

It is easy to see that \AF{y,t)\ = \F{y/2,t) - F{-y/2,t)\ < 2umax- Also \Au{y,t)\ < 
\y/t\ + \ AF{y — XQ,t)\ whenever one of the points y/2 or —y/2 is in the support of u. So we 
bound the first term on the right hand side of (49) by 

{i\y/t\ + \AF{y - xo)|)"X[,.^(,/2),r +{j//2)] (O)^(r* " *)) 
< 2"(|AF(y - xoTOiT* - t)) + 2"(|y/trx[,_(,/2),r +(j//2)] (O)^(T* - t)). (50) 

The first term on the right hand side of (^) can be bounded by averaging over xq first and 
using 



£^^\AF{y-xo,tr < (2nmax)"-^/'^^|AF(y-xo,t)| 

< (2nnrax)"'^^, 



using in the last inequality the fact that F is decreasing and bounded by 2tiniax- Substituting 
into ( |50D one can take the further averaging as for previous error bounds. By taking t large 
enough that L{t) > y and combining the various terms one arrives at the desired error bound. 
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6.5 The proof of the estimate (|35| ). 

The proof of this estimate is similar to that of (31). Noting that Au{y, t) = AF{y — xq, t) + 
{y/t) we may write 

j^e{u-Au{y,t)) = J^e{\AF{y-xo)\-{y/t)-\n\) 

dxo \AF{y-xo)\ 



< 



< 



2L {y/t) + \u\ 
1 2yumax 



2L {y/t)^\u\)' 
A similar estimate holds for Au*{y,t). Hence 



\Riu,y,t)\ < {i9iu-Au{y,t)) + 9iu-Au*iy,t))e{T*-t)) 

nax' 

L{t) y ftc^'/^ 



^ T ( Jll n (^maxg(T*-t)) 
L [y/t) + \u\) 



< 



L y +\u\\t 



6.6 Proof of the estimate (|H]). 

The proof of this estimate is similar to that of (^). The key change is to obtain a bound for 

J ^\F{y-xo,t)-F{y-xo,t + T)\. (51) 

The piecewise constant profile F{y, t) consists of a series of shocks which may travel forwards 
or backwards but move with a maximum speed Umax- The total height of the shocks is also 
bounded by Umax- So the integral (51) can be bounded by uj^g^-x^T/2L. The possibility of 
infinitely many shocks, or the merging of shocks between times t and t + r, does not affect 
this upper bound. 

6.7 The proof of the estimate (^. 

The construction of the two shock profile uses two parabolas that pass through the graph of 
the Brownian motion —q{x) = — J^^uo{z)dz at its point of maximum. Below is a lemma 
about the behavior of a Brownian path near its maximum. 

Lemma 1 Let (Bt : < t < 1) be a standard Brownian motion started at zero. Define 

M = sup Bt, S = mi{t : Bt = M}. 

te[o,i] 
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We consider the pieces of the path (Bt) either side of its maximum by defining 
Xt = M- Ss-t for t G [0, S], Xt = M- B^+t for t G [0, 1 - S]. 
Define the slopes of two lines that pass through the maximum and lie above the path by 
e = mi{Xt/t : < t < S}, G = mi{Xt/t : < t < 1 - S}. 

a) The triples {M, E, {Xt : t < S)) and (M — Bi, 1 — S, (Xj : t < 1 — S)) are identically 
distributed. 

b) The law of (Mj'E) is given by 

P{M G dm,T; £ da) = — — , exp(—m^ /2a)dmda. 

TT{a{l — a)YI^ 



c) Conditional on M £ dm, T, £ da the path {Xt : t < a) satisfies Xq = and solves the 
stochastic differential equation, driven by a Brownian motion (Wt), 

dXt = fit,Xt)dt + dWt, wi/iere/(t,x) = ^ + |^(exp(f^)-l)-i. (52) 

d) For (Xt) that solves (52) we have the estimate 

P(e <e)< C9{m + crm-i) + I(m < 9a). 



We delay the proof of this lemma until the end of this appendix and first use it to prove 
the estimate (0) on the tail P{T* > t) of the time T* at which the two shock profile is 
obtained. The construction of the two shock profile uses the function q{x) = uo{z)dz, 
its global minimum Q and the position xq + x* at which the minimum is attained. Two 
parabolas of the form tt{z) = {z — x)^/2t (and Tt{z) = P + {z — x)^/2t) are constructed to 
pass through the point (xq + x* , —Q). The slopes of the parabolas at the point xq + x* are 
(-2Q/t)i/2 (respectively (2(P - Q)/t)^/'^). Let T (respectively f) be the smallest time t at 
which the parabola tt (respectively vf) lies above the graph of —q{x). Then the two shock 
profile is attained for times t > T* = max{T, T}. 

To apply the lemma we must rescale to obtain a standard Brownian path of length 
one. Set Bt = —{2lJ)^^^'^q{xQ — I + 2lt) for t G [0, 1]. Then (Bt) is a standard Brownian 
motion and its maximum M takes the value —Q/{21J)^^^ . The construction of the parabola 
vr (respectively vf) show that if 6 < {-AlQ/tJ)^/^ then t < T (respectively if G < ((4/(P - 
Q)/tJ)^f'^ then t < T). Part a) of the lemma shows that both of these events have the same 
probability. So, applying part d) of the lemma, 

P(r* >t) < 2P(G < {-AlQ/tJ)^/^) 

= 2P(G < 2''i^M^'h-^/H^'\r^'^) 

< Ct-^lH^/^J-^/^E{M^/'^{M + T.M-^)) 

+2P{M < 25/4Mi/2rV2/3/4_^-i/45.^ 

< Ct-V2/3/4j-l/4 
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using Markov's inequality in the last inequality and the exact distribution of (M, S) in part 
b) of the lemma. This completes the proof of (^) and it remains to describe the proof of the 
lemma. 



Part a) of the lemma follows from the symmetry of the problem with respect to the time 
reversal i — > 1 — t. The distribution of (M, S) is well known and may be obtained for example 
by exploiting the reflection principle. Conditional on M G dm, S G cr the path {Xt) becomes 
a Brownian bridge, taking the value zero at time zero and the value m at time a, that is 
conditioned to never take negative values. The equation describing the evolution can then 



be obtained using an h-transform as in Rogers and Williams [30| section 4.23. 



We first sketch the idea for estimating P(Q < 0) = P{Xs < 9s for some s < a). The 
drift f{t,x) in equation (|52|) is approximately 1/x for small t and x. If this approximation 
were exact the process {Xt) would satisfy dX = X~'^dt + dW which is uniquely solved by the 
three dimensional Bessel process (the radius of a three dimensional Brownian motion). For 
a Bessel process one can make use of time inversion via the identity in distribution 

(Xt : i > 0) = {tX^,t : t > 0) 

and potential theory for three dimensional Brownian motion which gives 

P{Xs < for some s > 0\Xq = x) = mm{9x~^ , 1}. 

Then 

P{Xs < 9s for some s < a) = P{Xs < 9 for some s > 1/cr) 

= E{nun{9X-l^,l}) 

= E{mm{9a^/^X^\l}) 

= / (27r)-3/V2exp(-rV2)min{6'o-i/V-\l}dr 
Jo 

< C9a^l'^ 



where the penultimate equality follows from Brownian scaling and the final equality from 
a calculation using the density of the Gaussian variable X\. To exploit this idea we divide 
the interval [0,(j] into two parts, over the first of which the approximation f{t,x) ~ 1/x is 
sufficiently good. 

We first estimate P{Xs < 9s for some s < ct/2). Using the elementary inequalities (1 — 
z)/2z < (e^^ — < 1/22: for all z > one obtains the bounds x~^ — x{cr — t)~^ < f{t,x) < 
x~^ + 2ma~^ . Hence 

X^^ - lo'^Xt < f{t, Xt) < Xt^ + 2m{a - t)"^ for t < a/2. (53) 

So the solution of the equation 

dYt = Yt-^dt-2a-^Ytdt + dWt, Yq = 
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satisfies Yt < Xt for all t < t. To remove the unwanted —2a~^Ytdt in the drift of (Yt) we use 
a change of measure. Define a new probability measure Q by defining the Radon-Nicodym 
derivative M by 



M = 



dQ 
dP 



/2 



cr/2 



> 



exp(2cj^^ / 

exp(a-iy,% + 2a-2^'''V, 
exp(-3/2) 



Y.dWs - 2a' 

a/2 



a/2 



Y^ds) 



■ds- 3/2) 



The second equality here follows from Ito's formula. By Girsanov's theorem ( see |2^) the 
process (Yt) solves dY = Y~^dt + dW with respect to some Brownian motion (W) under 
Q, implying that (Yt) is a three dimensional Bessel process under Q. Writing Eq for the 
expectation under Q we have 

P(Xs < 6s for some s < a/2) < P{Ys < 9s for some s < a/2) 

= Eq(M-^I(Ys < 9s for some s < a/2)) 



< 
< 



e^/^Q{Ys < 0s for some s < a/2) 
COa^l^ 



(54) 



using the argument given above. 



It remains to estimate the probability P{Xs < 9s for some a/2 < s < a). We shall 
further condition on the value of X^/2- If ^a/2 ^ dr the evolution of {Xg : s £ [a/ 2, a]) is 
that of a Brownian bridge starting at r, ending at m and conditioned to take non-negative 
values. We write Qx for the law of a one-dimensional Brownian motion (Wt) started at x 
and we define Ha = inf{t : Wt < a}. Then, supposing r,q > 9a, we have 

PiXs < 9s for some s e [a/2, a]\X^/2 € dr) 
= 1 - P{Xs > 9s for ah s G [a/2, a][X^/2 G dr) 
< 1 - P{Xs > 9a for ah s £ [a/2, a]\X^/2 ^ dr) 
= 1 - QriHea > o/2\W„i2 G dm, > a/2) 

^ _ QrjHea > a/2,W„i2 £ dm) 
Qr{HQ> a/2,W,/2£dm) 

The reflection principle can be used to show that, for a < r,m. 



Qr{Ha > t,Wt £ dm) = (pt{m — r) — pt{m + r — 2a)) dm 

where pt{z) = (27rt)~^/^ exp(— z^/2t). Using this we rewrite the last expression as 

exp((m -|- r)'^/a) — exp((m + r — 29a)'^ /a) 
exp((m -|- r)^/a) — exp((m — r)^/a) 



(55) 
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(l-exp( ))-^4e{m + r -2r])exp{> i) \ ^ ) 



Thus 



a a 

for some G [0, Oct] by the mean value theorem 

— AuxT 

< (l-exp( ))~^4e{m + r - 2rj) 

a 

< Cd(l + -^){m + r) (using (l-e-^)-i < (7(1 + z-i)) 



PiXs < 9s for some s £ [t, a]\Xr G dr) 
< Ce{m + r + ar'^ + am'^) + I{r <ea)+I{m<ea). (56) 

We now undo the conditioning on £ dr. Using the upper bound in (53) and Ito's 

formula one obtains dXf < (3 + Ama~^ Xt)dt + 2XtdWt. Taking expectations one has 

E{Xf) < 3t + 4mCT~^ t E{Xs)ds 

Jo 

< (3 + mV-^)t + 40--^ / E{X'^)ds. 

Jo 

Applying Gronwall's inequality shows that E{X„/2) < {E{Xj^^)y/'^ < C(a^l'^ + m). By 
Markov's inequality P(X^ < Oa) < eaE{X~^). Using the comparison with a Bessel process as 
before we have E{X~^) < €^^'^Eq{Y~J:^ < Ca^^^'^. Using these bounds in (^) and combining 
with (^) leads to the estimate in part d) of the lemma. 
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